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Abstract. We study the problem of adsorption of self- interacting linear polymers 
situated in fractal containers that belong to the three-dimensional (3d) Sierpinski 
gasket (SG) family of fractals. Each member of the 3d SG fractal family has a fractal 
impenetrable 2d adsorbing surface (which is, in fact, 2d SG fractal) and can be labelled 
by an integer b (2 < b < oo). By applying the exact and Monte Carlo renormalization 
group (MCRG) method, we calculate the critical exponents i' (associated with the 
mean squared end-to-end distance of polymers) and (p (associated with the number of 
adsorbed monomers), for a sequence of fractals with 2 < & < 4 (exactly) and 2 < 6 < 40 
(Monte Carlo). We find that both i' and 4> monotonically decrease with increasing b 
(that is, with increasing of the container fractal dimension df), and the interesting fact 
that both functions, i^(&) and </>(&), cross the estimated Euclidean values. Besides, we 
establish the phase diagrams, for fractals with 6 = 3 and 6 = 4, which reveal existence 
of six different phases that merge together at a multi-critical point, whose nature 
depends on the value of the monomer energy in the layer adjacent to the adsorbing 
surface. 
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1. Introduction 

Statistics of a polymer chain in various types of solvents near an impenetrable wall 
(boundary) with short-range attractive forces has been extensively studied because of 
its practical importance [Q], and as a challenging problem within the modern theory 
of critical phenomena The most frequently applied model for a polymer chain 
has been the self-avoiding (SAW) random walk model (that is, the walk without 
self-intersections), so that steps of the walk have been identified with monomers that 
comprise the polymer, while the solvent surrounding has been represented by a lattice. 
These problems have been mostly studied for models situated on two-dimensional 
(2d) Euclidean lattices using various theoretical methods, such as the series expansion 
method, the renormalization group (RG) techniques, the mean-field approach, Monte 
Carlo simulations, and the conformal invariance method. In the last two decades these 
problems have been also studied in a case of fractal lattices embedded in the two- 
dimensional Euclidean space. On the other hand, in a more realistic three-dimensional 
case (for the Euclidean lattices, as well as for fractal lattices), a smaller number of results 
have been obtained. The study of fractal lattices has an advantage not only because 
their intrinsic self-similarity makes the problem more amenable to an exact approach, 
but also because these lattices as such may serve to model porous media. 

In this paper we report results of our study of a linear polymer situated in the 
three-dimensional (3d) fractal lattices that belong to the Sierpinski gasket (SG) family of 
fractals, assuming the interaction between two adjacent nonconsecutive monomers and, 
in addition, assuming adsorbing interaction with the walls of the fractal interior. This 
problem has been studied in the case of the 3d Euclidean lattices (see, for instance f^, 
and references quoted therein), but the number of obtained results is definitely smaller 
than in the corresponding 2d case. In the 3d case, the main endeavor has been manifested 
in attempts to establish phase diagrams 0, ^ in the interaction parameter space (which 
consists of the monomer-monomer interaction parameter and the adsorption energy 
parameter). In addition to the phase diagram, attempts have been made to calculate 
critical exponents that characterize various polymer configurations. The first work of 
this kind for a fractal was done by Bouchaud and Vannimenus 0, who applied the 
renormalization group (RG) technique for the 3d SG with scale parameter 6 = 2. Here 
we report results of our exact RG calculation for the 3d SG with 6 = 3 and 6 = 4, and 
our results for the sequence 2 < 6 < 40 obtained via the Monte Carlo renormalization 
group (MCRG) method. Therefore, this paper appears to reflect an effort to extend our 
previous studies performed in the case of the two-dimensional fractal lattices. Such an 
effort has been incited by the well-known fact that critical properties of a given model 
depend on dimension of the space in which the model is situated. 

This paper is organized as follows. In Sec. ^ we first describe the 3d SG fractals 
for general 6. Then, we present the framework of the RG method for studying the 
polymer adsorption problem on these fractals (taking into account the presence of the 
monomer-monomer interaction), in a way that should make the method transparent for 
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Figure 1. The first two steps (r = 1 and r = 2) of tlie self-similar construction of the 
3d SG fractals, in the case b = 4. 

exact calculations, as well as for the Monte Carlo calculations. In Sec. |] we elaborate 
on the phase diagrams obtained through the exact RG analysis for the 6 = 3 and 
6 = 4 SG fractals, and, in addition, we display our findings for the concomitant critical 
exponents u (associated with the mean squared end-to-end distances of polymers) and 
(associated with the number of adsorbed monomers). For 6 = 3, values of z/, for different 
phases, have been previously calculated 0, so that here we report the values for the 
corresponding crossover exponent 0, whereas in the 6 = 4 case we had to calculate values 
for the both exponents u and 0. In Sec. ^ we explain details of the MCRG calculations 
of the critical exponents, for arbitrary 6, and present their specific values up to 6 = 40. 
Summary of the obtained results and the relevant conclusions are given in Sec. ^. 

2. Renormalization group scheme 

Each member of the 3d SG family of fractals is labeled by the scale parameter 
6 = 2, 3, 4, . . . and can be constructed recursively starting with the pertinent generator 
G'^^\b) which is a tetrahedron of base 6, that contains 6(6+1) (6 + 2)/6 unit tetrahedrons 
(see Fig. [l|). The subsequent fractal stages are constructed self-similarly, by replacing 
each unit tetrahedron of the initial generator by a new generator. Thus, to obtain the 
rth-stage fractal lattice G^'^\b), which we shall call the rth order generator, the recursive 
process has to be repeated (r — 1) times, so that the complete fractal is obtained in the 
limit r oo. Fractal dimension dj of 3d SG fractal is equal to 

rff = ln[(6 + 2)(6 + l)6/6]/ln6. (2.1) 

We assume here that one of the four boundaries of the SG fractal is impenetrable 
attractive surface (wall), which is itself a 2d SG fractal with the fractal dimension 

rff = ln[6(6 + l)/2]/ln6. (2.2) 

In order to describe both the effect of monomer-monomer interaction and the effect 
of attractive (adsorbing) surface, one should introduce the three Boltzmann factors: 
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Figure 2. The fractal structure of the & = 2 3d SG fractal at the first stage of 
construction, with an example of the SAW path. The shaded area at the basis of 
the thetrahedron represents the adsorption wall. The steps on the adsorbing wall 
and in the adjacent layer are weighted by the factors w = e"*^™/"^ and t = e"*^'/"^, 
respectively. Here ew is the energy of a monomer lying on the adsorbing wall (e^, < 0), 
and > is the energy of a monomer that appears in the layer adjacent to the 
wall. The Boltzmann factor v = e.~'^^l^ corresponds to the energy of interaction 
El, < between two nonconsecutive neighboring monomers. The depicted SAW path 
represent one term in equation (2.5) for r = 1 with A^ = 5, P = 4, M = l and K = 2. 



V = e"'^"/'^^-^, w = e^^""/'^^-^, and t = Q-^t/ksT ^ where is the energy corresponding to 
interaction between two nonconsecutive neighboring monomers, is the energy of a 
monomer lying on the adsorbing surface, and St is the energy of a monomer in the layer 
adjacent to the surface (see Fig. 2). If we assign the weight x to a single step of the SAW 
walker, then the weight of a walk having N steps, with P nearest neighbor contacts, M 
steps on the surface, and K steps in the layer adjacent to the surface, is x^v^w^^t^ . An 
arbitrary SAW configuration can be described, following by using the five restricted 
generating functions (see Fig. 3). For G^'^\b), the generating functions, in terms of the 
interaction parameters, have the form 

A W {x,v) = Y^ A^"-^ {N, P)x^v'', (2.3) 

N,P 

5 W {x,v) = J2 ^^"^ (N, (2.4) 

N,P 

A^{\x,v,w,t)= J2 A'(\N,P,M,K)x''v''w^t'', (2.5) 

N,P,M,K 

A^;\x,v,w,t)= J2 At\N,P,M,K)x''v''w^t'', (2.6) 

N,P,M,K 

B[''\x,v,w,t)= 13^{\n,P,M,K)x''v''w^H'', (2.7) 

N,P,M,K 

where the coefficients have the following meanings: 

• A'^^^ {N, P) is the number of A^-step SAWs, lying completely in the bulk, with P 
nearest neighbor contacts, and entering G^'^\h) through one corner vertex, and 
leaving it via a second corner vertex, 

• B^'^\N, P) is the number of A^-step SAWs, traversing the G^'^\b) twice, in the bulk, 
with P nearest neighbor contacts, 

• A'{\n,P,M,K) {AP{N,P,M,K)) is the number of A^-step SAWs entering the 
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Figure 3. Schematic representation of the five restricted generating functions used 
in describing all possible polymer configuration within the rth stage 3d SG fractal 

(r) 

structure. Thus, for example, A\ represents the SAW paths that start at one 
tetrahedron vertex that lies on the adsorption wall, and exit at the other vertex that 
also lies on the adsorption wall. The interior details of the r-th order fractal structure 
are not shown (they are manifested by the wiggles of the SAW paths). 



through a corner vertex lying on the adsorbing surface, and leaving it via a 
second corner vertex on the surface (in the bulk, in the case of A2'\n, P, M, K)), 
with P nearest neighbor contacts, M steps in the surface, and K steps in the layer 
adjacent to the surface (see Fig. 2), and, finally, 

• B^\n, P, M, K) is the number of A^-step SAWs going twice through the G^'^\b), 
with P nearest neighbor contacts, M steps in the surface, and K steps in the layer 
adjacent to the surface. 

These generating functions (depicted in Fig. ^ are parameters in the renormalization 
group (RG) approach, and for any b >2, RG equations have the form 

Air+i) = J2 a{NA,NB)A^^B^^, (2.8) 

Na,Nb 

= ^ b{NA,NB)A^^B^^, (2.9) 
A^'^ = E ai(Ar^,Ar^,Ar^„iV^„iVBjA^^5^M^Mf^^S^\ (2.10) 

Na,Nb,Na^,Na2,Nb-^ 

At'^ = E a,{NA,NB,NA,,NA,,NBM'''B''-Af^^A^^^B^^\ (2.11) 

Na,Nb,Na^,Na2,Nbi 

S^'^ = E b,{NA,NB,NA,,NA,„NBM'''B''-Af'^A^^^B^^\ (2.12) 

Na,Nb,Na-,,Na2,Nb^ 

where we have omitted the superscript (r) on the right-hand side of the above 
relations. The self-similarity of the fractals under study implies that numbers 
a{NA,NB), b{NA,NB), ai{NA, Nb, Na„ Na„ NbJ, a2{NA, Nb, Na„ Na„ NbJ, and 
bi{NA, Nb, Nai, Na2, NbJ, of the corresponding SAW configurations within the 
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structure do not depend on r. Starting with the initial conditions 

A^°\x, v) =x + 2x^v + 2x\^ 

A^i^x, V, w, t) =wx+ (w^ + t'^)x'^v + 2wt^x^v^, 



A^2\x, V, w, t) =tx + 2twx^v + 2tw^x-^V^, 

b[^\x, V, w, t) = wtxh^, (2.13) 

which correspond to the elementary tetrahedron G^^\b), one can iterate RG relations 
(p.8|)-( prT2D for various values of interactions v, w and t, and explore the phase 
diagram. This approach (which implies that interactions are restricted to sites within 
the elementary tetrahedron G^^^ and, moreover, that SAW exits G^**-* whenever it reaches 
its corner vertex 0]) was applied in ^ for 3d 6 = 2 SG fractal. Here we present an 
analogous type of analysis for the larger b cases. 

The average number of monomers in contact with the adsorption wall, for SAW 
spanning a G^^^ (6) can be expressed in terms of the partial derivatives of the generating 
functions Ai''^ and Ag''^ 

J:n,p,m,k M {At^ {N, P, M, K) + ) {N, P, M, K)) x'^vVH'' 



w 



OA'; 



A^;^ + aP 



A^{^ + aP 



+ 



w 



4''^ + A 



(r) 

w 



dw dw 

whereas the total average number of monomers can be expressed in the form 
En,p,m,k N {At\N, P, M, K) + At\N, P, M, K)) x^vVt^ 



(2.14) 



X 



4'^ + A^^^ 



dx 



At^ + A^^ 



+ 



X 



dx 



A^^^ + At 



- (4:1 + 41) • (2-15) 



From the RG equations ( p.8| )-( prT^ ) one can obtain recursion relations for the derivatives 
of the generating functions in the following matrix form 



(2.16) 



where matrices Rs and R are comprised of partial derivatives of generating functions 
A\ B', A[, A'2 and corresponding to SAWs spanning the generator G^''"''^^ (6) , in 
respect to generating functions A, B, Ai, A2 and Bi, corresponding to G^^'\b): 
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Starting with the initial conditions for the derivatives 



dx 
dx 



l + 2xv + Qx^v^ 



dx 



2xv^ 



(0) 



dx 
1? 



Mi") 



w + 2{w'^ + t^)xv + Qwt^x'^v^ 

— = X + 
dBf^ 



2wtxv 



dw 



dw 

2tvx'^ + Atwx^v^ 



a4f 

dx 



X + 2wx^v + 2t'^x^v^ 



t + itwxv + Qtw^x'^v^ , 



dw 



txh^ 



(2.17) 



and iterating recursion equations (|2.16| ), it is possible to establish the relation between 
the average number of adsorbed monomers (M('')) (gJ§ and the average length of the 
polymer chain {N^^^) (|2.15| ) in the limit r oo, for various values of the interaction 
parameters v,w, and t. 



3. Exact approach: phase diagram and critical exponents 

To solve exactly the adsorption problem of a self-interacting SAW for arbitrary member 
(for any 6 > 2) of the 3d SG family, it is necessary to find all coefficients that appear in 
the set ( |2.8D -( P?T^ ) of the RG equations. The 6 = 2 case has been completely analyzed 



in P] , while in the 6 = 3 and 6 = 4 cases only the bulk RG equations (p.8|) and (p.9|) have 
been studied, in and [§] respectively. In this paper we make a complete analysis of 
the 6 = 3 and 6 = 4 cases, that is, including the adsorption RG equations (|2.1CI| )- (|2.12|) . 

3.1. Theb = 3 SG fractal 

The RG equations (|2.(j|) and (|2.9|) for the bulk parameters A and B, found in H] for the 



case of the 6 = 3 SG fractal, have the form 



^(r+i) = ^3 ^ g^4 ^ ^g^5 ^ 34^6 ^ 7g^7 ^ ^-^2A^ + 112A^ + 64^4^° + 8A^B + 
36A^B + UOA^B + 292/l^S + 424^^5 + 332^^5 + 12A^B^ + 12A^B^ + 
118^^5^ + 380^^5^ + 806^^5^ + GGiA^B^ + 72A^B^ + 352^^5^ + 
704^^5^ + 1728^^5^ + UAA^B^ + 1568^^5^ + 848^^5^ + 

264A^S^ + 3192A^S^ + 320/l=^E^ (3.1) 
^(r+i) = ^6 ^ ^2A' ^ 40^8 ^ gg^9 ^ 32^10 ^ 28A^B + 88A^5 + 224^^5 + 

imA^B + 40^*^5^ + A9QA'^B^ + 596^^5^ + 176^^5^ + 768^*^5^ + 
1056^^5^ + 88A^B^ + 2UA^B^ + 2bUA^B^ + 1152A^S^ + 

1888^^5^ + 5808/1^5^ + 1936^^5^ + 4308^^5^ . (3.2) 

Here we give a summary of the analysis of the above set of equations. For small 
values of the monomer-monomer interaction {v < vq), the extended SAW phase fixed 
point {Ae,Be) = (0.34196, 0.02395) is reached through the RG transformations. 
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Linearization of the RG equations in the vicinity of this fixed point gives only one 
relevant eigenvalue = 5.36201. The mean-squared end-to-end distance {R%) of the 
A^-step polymer chain, in general case, behaves asymptotically (for >> 1) as 

{R%) ~ N'", (3.3) 
where the critical exponent u is given by 

which for 6 = 3 gives ue = 0.6542. Starting with v = vq, the RG equations lead to the 
fixed point {Ag, Be) = (0.20717, 0.43075), for which both eigenvalues = 8.72308 and 
A^ = 2.45012 are relevant. At this fixed point critical exponent z/^ = In 3/ In A^ = 0.5072 
is smaller than ue, which is the manifestation of the so-called collapse transition. At 
the collapse transition, the free energy per site / behaves as 

/ ~ |t; - vgf~°' , (3.5) 

where the critical exponent a in general case is given by 

a = 2-^, (3.6) 

and in this specific case it is negative, a = —0.4170. Depending on the value of 
the one-step weight (fugacity) x, for strong monomer-monomer interactions {v > vg) 
RG equations ( |3.1|) and ( p.2| ) bring about the trivial fixed point {A,B)* = (0,0), for 
X < x*{v), or {A, B)* = (oo, oo) for x > x*{v), whereas for x precisely equal to x*{v), the 
fixed point {Ac, Bq) = (0, oo) is reached. Analyzing the RG equations in the vicinity of 
{Ag, Bg), by keeping only the dominant terms in the right-hand side of ( |3.1| ) and ( |3.2| ), 
one can find = 0.48195. Since the fractal dimension dY^^ = I/^'g = 2.07491 of the 
SAW in this case is larger than the fractal dimension of the extended SAW (for v < vg), 
but still less than the fractal dimension of the underlying lattice df = 2.09590 one 
can conclude that fixed point (0, cxd) corresponds to the 'semi-compact' phase This 
finding is in contrast with the results obtained for polymers on homogeneous lattices 
and onthe6 = 23dSG fractal lattice 0, where one finds d^^^ = df. 

In order to establish the specific form of the complete set of the exact RG equations 
(p.8|)-( p?I^ ), required for the study of the adsorption problem on the 6 = 3 SG fractal, 
we have enumerated the requisite SAW configurations, achieving thereby the pertinent 
RG coefficients. This procedure is rather complex, as well as the corresponding set of 
coefficients, and for this reason we give them in the Appendix A. The physical picture 
that follows from these RG equations, for various values of the interaction parameters 
{v,t and w), is unusually rich and we present it in what follows. 

Numerical study of the adsorption RG equations (given in the Appendix A) 
shows that for w > 1 (attractive surface) an unbinding transition appears at a finite 
temperature only if t < 1 (repulsive interaction in the layer adjacent to the surface). 
The nature of this transition depends on the value of the monomer-monomer interaction 
V. This situation is analogous to the previously studied cases of the two-dimensional 
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and three-dimensional 6 = 2 Sierpinski gasket fractals P, but the new physical features 
do appear. 

3.1.1. Extended SAW phase For the chosen initial conditions ( |2.13| ) we found the 
critical value of monomer-monomer interactions ve = 2.446161 (which is different from 
the value found in because of the slightly different initial conditions, but which 
does not affect the overall critical behavior). For weak monomer-monomer interactions 
V < vg, and the corresponding critical fugacity x = x*{v), RG equations ( p.l|) and (|3.2| ) 
for the bulk parameters A and B iterate towards the fixed point {Ae, Be)- 

Behavior of surface RG parameters {Ai \ A'2 \ B^"^) depends primarily on the 
corresponding interaction parameter w. Thus, for weak interactions w < Wc{t,v), 
and X = x*{v), the parameters [A^p, A^2\ ^1^^) tend to (0,0,0), when r —>■ oo, which 
indicates that polymer, being in the extended coil phase, stays away from the attractive 
surface. We shall refer to this state as the desorbed extended (DE) phase, determined 
by the RG fixed point 

(A B, A,, A2, B,)* = {Ae, Be, 0, 0, 0) . (3.7) 

This behavior changes abruptly at w = Wc {t, v), where A^^A''^^ Ae and b[''^ Be- 
The new behavior is described by the new fixed point 

{A, B, Ai, A2, BiY = {Ae, Be, Ae, Ae, Be) , (3.8) 

which is the symmetric special fixed point, that corresponds to the unbinding transition 
of the SAW. If w is increased beyond Wc{t,v), and for the fugacity equal to x = x*{v), 
bulk parameters still approach {Ae,Be), but the surface RG parameters diverge, 
implying that critical fugacity Xc is smaller than x*{v). A thorough analysis shows 
that Xc depends on the values of v, w, and t, whereas RG parameters flow towards the 
new fixed point 

{A, B, A,, A2, B,)* = (0, 0, Af, 0, 0) , (3.9) 

where A"^ = 0.551147 is the fixed point for the 6 = 3 two-dimensional SG fractal (see 
0]), meaning that for strongly attractive surface polymer remains adsorbed (this is the 
adsorbed SAW phase). 

At the symmetric special fixed point ( |3.8|) linearized RG equations have two relevant 
eigenvalues: the larger = 5.36201, already found for the bulk RG equations, and the 
smaller Af = 3.32923. The matrix of the linearized RG equations is equal to the matrix 
R in ( |2.16| ), which means that 

(ArM)~(Af)', for r>l, (3.10) 

whereas the largest eigenvalue of the matrix R5 in ( 2.16 ) is equal to Af , meaning that 



(M^^)) ~ (Af)' , for r>l. (3.11) 
Consequently, for r — > 00, it follows 

(M(")) ~ (iV^^y , (3.12) 
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where the crossover exponent (f) is equal to 

In Af , , 

0£; = ^ = 0.7162. (3.13) 

This completes our analysis of the extended polymer phase. 

3.1.2. 9-line (v = vg) Here we organize our discussion of various phases that appear 
along the 6'-line in the phase space {v, w), for v = vq = 2.446161 and for various values 
of w. 

For X = x*{vg) = 0.109683 and w < wg{t), RG iterations lead to the fixed point 

(A B, A,, ^2, B,y = {Ae, Be, 0, 0, 0) , (3.14) 

which means that for small values of w polymer remains dcsorbed, in the solution, in 
the form of the ^^-chain. On the other hand, for w > We{t) critical fugacity Xc{v0,'w,t) 
is less then a;*(fg), and the relevant fixed point is again (0,0,74^,0,0) which describes 
an adsorbed two-dimensional SAW (as in the case v < vq). 

At the critical value w — WQ{t), and x still equal to x*{vq), RG parameters flow 
towards the new fixed point 

(A, B, Ai, A, BiY = {Ae, Be, A^e, A^e, Bw) , (3.15) 

whose coordinates Aie,A2e, and Bie depends on the value of t, and accordingly there 
are five possible situations that should be analyzed. 

First, for < t < t* = 0.1553901, the new fixed point is 

{A, B, A, A2, B,y = {Ae, Be, Af, 0, 0) , (3.16) 

that controls the coexistence between the ^-chain in the bulk and the adsorbed two- 
dimensional SAW. 

Second, we find the fixed point 

(A, B, Ai, A2, BiY = (Ae, Be, 0.5234, 0.1637, 0.0305) , (3.17) 

for the critical value of the interaction t — t^, with four eigenvalues larger than one 
— the two bulk values and A^, and the two surface eigenvalues A^^ = 3.96291, and 
A52 = 1.28383, which brings about 

05 = ^^ = 0.6357. (3.18) 

Third, in the interval between the two critical values of the interaction parameter 
t,tl<t<t*2^ 0.6573781, the new fixed point is 

(A, B, Ai, A2, Bi)* = (Ae, Be, 0.19265, 0.44713, 0.132446) , (3.19) 

for which, in addition to the relevant bulk eigenvalues, there is only one surface relevant 
eigenvalue A5 = 3.663797, and the crossover exponent is equal to 

= rrf = 0-5995 . (3.20) 
In A?, 
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Fourth, for the second critical the RG parameters flow towards the symmetric 
fixed point 

(A B, A,, A2, B^y = {Ae, Be, A, A, Be) , (3.21) 

where, as in the case of the first critical value t\, one finds two relevant surface eigenvalues 
-^si ~ 5.368208 and = 1.718244, wherefrom one obtains the crossover exponent 

<^^ = 7^ = 0-7758. (3.22) 
In \l 

Fifth, for ^2 < < 1 the following fixed point 

(A 5, A, A, BiY = {Ae, Be, 0, 0, cx)) (3.23) 

is reached. By keeping the dominant terms on the right-hand side of the RG equations 
( |A.1| - |A.3| ), for the surface parameters Ai,A2, and Bi (see Appendix A), the RG 
equations attain the approximate form 

A[ ^ {128AI + 124A^ + 292A^ + 264^^5^ + 944^^5^ + 320AgB^)AlB^ = CiAlB^ , 

B[ ^ (472A^ + U52AlBg + 1452^^5^ + 4308A5^)Ai5f = C2A1BI , 

A'2 ^ C3A2 , (3.24) 

where ci ~ 6.26, C2 ~ 89.32, and C3 ~ 0.127. Introducing the new variable yi = BiA^i^ 
we obtain the tractable form of the approximate RG equations 

A = ciy^Ai , y[ = cl'^C2y\ , A = C3A2 , (3.25) 

which have the fixed point (0,0,?/^ = Ci ^^^°C2 ^''^). Linearizing these RG equations 
in the vicinity of the fixed point, one relevant eigenvalue A5 = 6 is found, with the 
corresponding crossover exponent 

In A? 



InAf, 



0.8272. (3.26) 



3.1.3. Semi-compact regime (v > ve) For values of the monomer-monomer interaction 
parameter v larger than ve, and for the fugacity x equal to the bulk critical value 
x*{v), the bulk RG parameters {A,B) fiow towards (0, 00), while critical behavior of 
the surface RG parameters {Ai, ^2,-81) depends on the values of both surface interaction 
parameters w and t. In particular, for any value of t between and 1, there is a critical 
value Wc{t, v) such that {Ai, A2, Bi) fiows towards (0, 0, 0), for w < wdt, v), whereas for 
w precisely equal to Wc(t,v) one observes 

(A, B, A,, A,, BO ^ j ~- -^l^- °f • < "» + ' (3.27) 

[ (0, cxo, 0, 0, cxd) , V > Ve + e 

for t <t* ^ 0.156, and 

{A, B, Ai, A2, Bi) (0, 00, 0, 0, 00) , (3.28) 

for t > t* and all v > ve- 

The fixed points (0, 00, 0, 0, 0), {0, 00, Ag^^^r,0,0), and (0, 00, 0, 0, 00), correspond 
to the desorbed semi-compact polymer chain (globule), to the coexistence between 
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the globule in the bulk and the adsorbed 2d polymer, and to the adsorbed globule, 
respectively. Finally, for values of the parameter w larger than the critical value Wc{t, v), 
the critical fugacity Xc = x*{v, w, t) is smaller than its bulk critical value x*{v), and RG 
parameters iterate towards the fixed point (0,0,741^^,0,0), which corresponds to the 
fully adsorbed polymer. 

The adsorbed globule phase can be analyzed by keeping the dominant terms on 
the right-hand side of the RG equations (|3.1| ), ( ^72|) , and (|A.1| )- (|A.3| ), in the vicinity of 



the fixed point (0, oo, 0, 0, oo). Accordingly, one can obtain the following approximate 
equations: 

A' ^ 320A^B^ , B' ^ 4308^^5® , A[ ^ 320AAIB^B^ , 

A'^ ^ AAAAfBBf , B[ ^ A308AAiB^Bl . (3.29) 

Introducing new variables 

y = AB\ yi = Ar/A, ^2 = ^7^2, ys = BiA\ (3.30) 

where 



73-5 I9 + V/73 

Q= 77; , (3-31) 



12 ' ^ 12 
equations (p.29|) transform into more tractable form 



A' = 320A^~Y , y' = 320^4308 /+^" , (3.32) 
y'i = ^^ 2/3 = 320^4308 yiy^'^'yl , (3.33) 

y 



y2 = -^'-rj- (3-34) 



J 320"? 1/6^-^ 
44 yfyl 

The new equations ( ^^321) have fixed point A* = 0, y* = (4308 x 320^)"^/^^+^^), with one 



relevant bulk eigenvalue = 8+Qz = (ll-|-\/73)/2 (see ^j). Inserting y* into equations 
( p.33| ) one finds that corresponding equations for the fixed point {yl,y^) are mutually 
linearly dependent. More precisely, the only fact that springs from these equations is 
the relation 

y* = {yiy^'yl- (3.35) 

For various values of t, large enough v and corresponding bulk critical fugacity x*{v), 
and the critical value Wcit,v), point {yi.yz) tends to different fixed points {yl,y^), but 
the above relation stays satisfied. Obviously, knowing {yl, y^), from equation (|3.34|) one 
can calculate the fixed value 2/2, which is also different for various t and v, but on the 
other hand is in excellent agreement with the value obtained via explicit iteration of the 
expression A'^/A2. 

Linearizing equations (p.33| ), ( |3.34| ) around the corresponding fixed point {yl, y2,yl) 

to 

In 6 



the second relevant eigenvalue Af = 6 is found, so that the crossover exponent is equal 
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This result is in agreement with the value estimated via direct numerical analysis, 
performed using equations ( |2.14| ) and ( p.l5| ), as explained in Section One should 
also observe that obtained value of 0g is slightly larger than dY/dY = 0.7781, which 
has been predicted in |Q for SAW in the compact phase. Of course, due to the topological 
frustration, SAW on 3d 6 = 3 Sierpinski fractal lattice can not, even for large monomer — 
monomer interaction, have a compact configuration (see ||^). Instead, it is in the semi- 
compact phase, in which its fractal dimension is less than the fractal dimension of the 
lattice (although larger than fractal dimension of SAW for v < vq). 



3.2. 3dh = A SG fractal 

As the scaling parameter h increases the number of possible polymer configurations on 
the 3d SG fractal lattices quickly grows. The RG equations, for 6 = 4, for the bulk RG 
parameters A and B have been found and analyzed by Maricic and Elezovic-Hadzic (see 
and Appendix B). The conclusion has been reached that qualitatively the physical 
picture of the polymer behavior in the bulk, for 6 = 4, is similar to the cases 6 = 2 and 
6 = 3. There are three possible phases in which polymer can reside, which we review in 
the following paragraph. 

For small values of the monomer-monomer interaction v < Vg = 2.33187 polymer is 
in the swollen state. The fixed point (A, 5)* = {Ae,Be) = (0.2899,0.0122) is reached 
for any v < Vg and for the corresponding critical fugacity x = x*{v). Linearizing 
RG equations around this fixed point one gets single relevant eigenvalue = 8.6924, 
and the concomitant critical exponent u is equal to ue = ln6/lnAf = 0.6410. On 
the other hand, the low-temperature fixed point (A, B)* = (Aq, Bq) = (0, 22~^/^) is 
reached when the RG iteration starts with v > vq. In this case, there is one relevant 
eigenvalue = 16, and the end-to-end critical exponent is = 1/2- The fractal 
dimension of the chain d^^^ = l/z/^ = 2 is less than the fractal dimension of the lattice 
df = In 20/ In 4 ^ 2.16. This means that polymer is in the semi-compact phase for 
strong monomer-monomer interactions, as in the 6 = 3 case, and in contrast with 
the cases of polymer on the 6 = 2 fractal and on the Euclidean lattices. Finally, 
when V = vg, the tricritical fixed point {A, B)* = {Ae, Be) = (0.1929,0.3388), that 
corresponds to the collapse transition, is reached. In this case there are two relevant 
eigenvalues, A^ = 15.4230 and A^ = 5.5357. Consequently, the critical exponent z/ is 
equal to ug = In 6/ In A^ = 0.5067. In addition, we have obtained the specific heat critical 
exponent a = 0.4012 > which reveals the singular behavior, which was observed in 
the case 6 = 2 but not in the case 6 = 3. 

Recursion relations for the surface RG parameters Ai, A2, and Bi are cumbersome 
(each of the corresponding equations has more than 3000 terms) and we are not 
going to quote them here (but they are available upon the request to the authors). 
Detailed numerical analysis of these RG equations shows that, depending on the value 
of interaction parameters v, w, and t, the following phases and the corresponding fixed 
points are accessible: 
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• extended desorbed phase, {Ae, Be, 0, 0, 0), for v < vg and w < Wc(t, v), 

• semi-compact desorbed phase, (0, Bq, 0, 0, 0), for v > vg and w < Wc{t, v), 

• desorbed 6'-chain, {Ag, Bg, 0, 0, 0), for v = vg and w < wg(t), 

• fully adsorbed chain, (0, 0, A^^^y, 0, 0), for w > Wc{t,v), where ^g^iy = 0.5063 is 
the fixed point for the 2d 6 = 4 SG fractal (see 0), 

• coexistence of the globule in the bulk and the 2d adsorbed polymer chain, 
(0,-8^5^^5 0,0), for Vg < V < Vg + e (where e is a small positive number whose 
specific value depends on t) and w = Wc(t,v), 

• surface attached extended chain, special symmetric fixed point, {Ae, Be, Ae, Ae, Be) 
for V < Vg and w = Wc(t, v), with the relevant eigenvalues and Af = 4.4533 that 
give the crossover exponent (j)E = InAf/lnAf = 0.6907, 

• surface attached globule, (0, Bq, 0, 0, Bq), ioi v > vg + e and w = Wc(t, v), with the 
relevant eigenvalues A^ and Af = 9 and (pc = 0.7924, 

• surface attached 6'-chain, multi-critical point {Ag, Bg, Aig, A2g, Big), for v = vg and 
w = wg(t), with 

r (0.2275, 0.3659, 0.0433), t < t* 

{Aie,A2g,Big) = I {Ag,Ag,Bg), t = t* = 0.9577 (3.37) 

[ (0,0,0.3643), t*<t<l 
for which one finds 

r 4.9036, t<t* r 0.5811, t<t* 

A^ = I 8.4078, t = t* and 0e = J 0.7782, t = t* (3.38) 

[9, t* <t<l [ 0.8031, t* <t<l 

In the last part of this section we summarize, in Table [^, the numerical results 
obtained via the exact RG approach, for the 3d SG fractals b = 2,3, 4, and provide the 
relevant discussion of the pertinent phase diagrams. One should notice that we first give 
the fixed point values of the bulk parameters A* and B* and the accompanying end-to- 
end critical exponent u. For each b there are three fixed points, which correspond to the 
three bulk phases — the extended polymer phase, the ^-chain phase, and the globular 
(collapsed) phase. Furthermore, for a given b, one may observe that A* decreases, while 
B* increases, with increasing of the monomer-monomer interaction v, and consequently 
u decreases. On the other hand, when b increases u decreases for the extended phase 
and the 6'-chain phase, whereas in the case of the globular phase u does not display a 
monotonic behavior. Besides, in the globular phase, for 6 = 3 and 6 = 4 the exponent 
u is larger than the reciprocal of the fractal dimension of the underlying lattices, which 
implies that the globular phase is not compact. This is in contrast with the 6 = 2 case 
(as well as, in contrast with the Euclidean lattices), where u = l/dj. 

In Table |l| we also give the fixed point values of the surface RG parameters 
74^,^42, Bl and the corresponding values of the crossover exponent 0, for the unbinding 
transitions from adsorbed polymer phase to various desorbed phases. The relevant phase 
diagram is given in Fig. 4, for the 6 = 3 SG fractal (a similar phase diagram for 6 = 2 
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Table 1. Fixed points and critical exponents v and </> for different self-avoiding walk 
phases on 3d Sierpinski fractals, obtained via exact renormalization group approach. 
Values for the bulk fixed point {A* , B*) and the critical exponent for b = 2 and 
3 fractals were obtained in Q and |^ respectively, whereas the surface bulk point 
{A^ , A2 , Bl ) and the crossover exponent 4> were previously known for the b — 2 fractal 
only §. 
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was obtained in 0, whereas in the course of this work we obtained a diagram, of the 
same type, for 6 = 4). In Fig. 4, the unbinding transitions are represented by the curve 
that separates the "adsorbed chain" region from the "desorbed chain" regions. On this 
curve hes the multi-critical point, so that the part of the curve for smaller values of 
the parameter v {v < ve) corresponds to the extended attached chain phase, while for 
larger values oi v {v > vq) the curve corresponds to the attached semi-compact phase. 
The fixed point that defines the extended attached chain phase is a symmetric special 
fixed point, that is, ^4^ = ^42 = A* and B\ = B*, with the corresponding crossover 
exponent which decreases with increasing of the scaling parameter b (see the top part of 
the Table |l]). Here we would like to mention that this behavior of the crossover exponent 
will be demonstrated for larger b as well, using the MCRG method (see the next section 
of this paper). 

Continuing our discussion guided by the results presented in Table |l], we focus now 
on the attached globular phase {v > vg). For this phase, the fixed point parameters 
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satisfy the relations Al = = A* = and = B*, where B* = Tl'^l"^ for 6 = 2 and 
6 = 4, while 5* = oo for 6 = 3. It is interesting to observe that the critical exponents 
V and are the same for the fractals with 6 = 2 and 6 = 4, while for 6 = 3 these two 
exponents are somewhat smaller. Besides, we should point out that only for 6 = 2 the 
crossover exponent is equal to the ratio (i^'^/ S"^^ which is related to the fact that only 
for 6 = 2 the globular phase is compact. Finally, by comparing values of the crossover 
exponent 0, we may observe that in the attached globular phase the number of adsorbed 
monomers is relatively larger than in the corresponding Euclidean phases (0 = 2/3). 

Results obtained for the multi-critical points for the fractals with 6 = 3 and 6 = 4 
(see the middle right part of Table I) appear to be dependent on the parameter t that 
describes energy of a monomer in the layer adjacent to the adsorbing boundary. This is 
manifested by several possible fixed points, in contrast to the case 6 = 2. Particularly, 
the position of the multi-critical point, as well as the shape of the critical line Wciv) 
(that separates the adsorbed phase from the desorbed phases), depend on the particular 
value of t (see Fig. For the case 6 = 3 there are four fixed points that correspond 
to the surface attached 6'-chain, two of which (the first and the third one) appear to be 
very repulsive (each with four eigenvalues larger than one) and can be approached for 
the critical values t\ = 0.1553 and = 0.6573 of the parameter t. The latter fixed point 
appears to be symmetric (A^ = A2 = A* = Aq and B^ = B* = Bg) and accordingly is 
the only point that describes the isotropic chain. The other two fixed points, the second 
and the fourth, are less repulsive (each of them has three eigenvalues larger than one), 
and can be reached for any value of the parameter t in the intervals < t < ^2 ^ind 
^2 < t < 1, respectively. 

Contrary to the expectation, that may be formed on the findings presented for the 
cases 6 = 2 and 6 = 3, that number of attainable multi-critical points will continue 
to increase with increasing 6, in the 6 = 4 case only three multi-critical points can be 
reached. Indeed, for 6 = 4 there is only one critical value t* = 0.9577 of the parameter 
t such that for < t < t* and for 1 > t > t*, two different nonsymmetric fixed points 
(one for each of the two intervals) with three eigenvalues larger than one are reached. 
On the other hand, highly repulsive symmetric point, with four eigenvalues larger than 
one, can be approached only for t = t*. Each one of the three multi-critical points 
describes its pertinent surface attached ^-chain which is manifested by the facts that 
Al> Al> A* = Ae and 5* < Bg for t < t*, A* = A* = and 5* < Bg for t > t% 
meaning that the corresponding two phases are anisotropic, whereas for the critical value 
t = t* the surface attached 6'-chain is isotropic. Finally, one can notice that crossover 
exponent increases with t (which was not the case for 6 = 3), and moreover that the 
particular values of for 6 = 3 and 6 = 4, in their relevant intervals 1 > t > ^2 
1 > t >t*, are rather close. 
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Figure 4. Phase diagram in the case of the 6 = 3 SG fractal for two different values 
of the interaction parameter t: (a) t = 0.1 and (b) t = 0.8. For t — 0.1, close to 
the multi-critical point on the critical line exists a small region (vg < v < vg + e) 
of the first order phase transitions (where the globule in the bulk and the adsorbed 
chain coexist) which ends at the point depicted by the small dash. Beyond the dash 
(v > Vg + e) the transition at the critical line is continuous, and corresponds to the 
attached semi-compact chain. 



4. Monte Ccirlo renormalization group calculation 

In this section we are going to apply the Monte Carlo renormalization group (MCRG) 
method to calculate the critical exponents u end for 3d SG fractals with 6 > 5. First, 
we shell present the MCRG calculation of the critical exponent u. In order to find 
the SAW critical exponent u, in the bulk phase, we should determine the nontrivial 
fixed points of the RG transformations (|2.8|) and (|2.9| ), and then we should solve the 
corresponding eigenvalue problem of the linearized RG transformations, that is, we 
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should solve the equation 

dA' * 

= 0, (4.1) 



' dA' _ \ \ dA;_ 
> dA -^^z dB 

dB' ( dB' _ \ 

dA V dB 



where the asterisk means that all derivatives should be taken at the corresponding fixed 
point and the superscript prime is used instead of the superscript (r + 1). Knowing 
the relevant eigenvalue X,y, we can determine the critical exponent u using the formula 
( p. 41 ). To learn a specific value of ly, for a given b, one should first find the coefficients 
a{NA, Nb), and b{NA, Nb) in the RG equations ( |2.8| ) and (|2.9|). As it was detailed in the 
previous sections, it has been possible to calculate the exact values of these coefficients 
only for 6 < 4. Thus, to get an entire sequence of values of u for 6 > 5, we are going 
to circumvent the problem of explicit determination of the exact coefficients in the RG 
equations by applying the MCRG technique. 

The MCRG method, allows direct calculation of derivatives that appear in the 
eigenvalue equation (|4.1|) . It starts by locating the bulk nontrivial fixed points, which 
requires, at the beginning, implementation of a MC simulation of the SAWs for a chosen 
initial set of values {Aq, Bq). In other words, we let the walker start his walking, at one 
fixed corner of the fractal generator and record the other corner, at which it leaves the 
generator, together with recording the total numbers A^^ and A^^, of crossings of the A, 
or B, type through the elementary tetrahedron. The SAW walker crosses an elementary 
tetrahedron in the A (or B) way (see Fig. ^ with the weight (probability) Aq, and Bq, 
respectively. We repeat this MC simulation L times, for the same set (^40,-80)- Thus 
we find how many times the walker has passed through the generator in the A (or B) 
way and by dividing the corresponding numbers by L we get the values of the functions 
(PD and (^), denoted here by A'{Ao,Bo) and B'{Ao,Bq). 

In this way we get the value of the sums ( p.8| ) and ( p.9| ) without specifying the 
coefficients a(A'^, A'^), and b{NA, Nb). Then, the subsequent values An and Bn {n > 1), 
at which the MC simulation should be performed, can be found by using the generalized 



"homing" procedure [TI], [1^, which can be terminated at the stage when the 
differences An — An-i and Bn — Bn-i become less than the statistical uncertainties 
associated with An-.i and respectively. Consequently, fixed point {A*,B*) can be 

identified with the last value {An, Bn) found in this procedure. 

Having learnt the fixed point, we need to solve the eigenvalue equation (|4.1| ) in 
order to find the critical exponent u via the formula ( p.4|) . Thus, we should find the 
partial derivatives dY'/dX (where X,Y & {A,B}) at the fixed point. For instance, 
starting with ( |2.8|) (in the notation in which the superscript prime is used instead of the 
superscript (r + 1)) and by differentiating it with respect to A we get 

^= E NAa{NA,NB){A)''^'\B)''- . (4.2) 

Treating now A' as the grand canonical partition function, for the ensemble of all possible 
SAWs that start at one fixed corner of the fractal generator and leave it at the other 
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one, we can write the corresponding ensemble average 

{Na{A,B))a' = ^ E NAa{NA,NB){Af-{Bf-, (4.3) 

which can be directly measured in a MC simulation. Finally, comparing the last two 
equations, we can express one of the requisite derivatives in terms of the measurable 
quantity 

BA' A' 

^ = ^(iVA(A,i?))^,. (4.4) 

In a similar way one can find additional three derivatives, so that we can write the 
general formula 

dY' Y' 



ax A--<^->'- (*-^) 

where X and Y stand for any pair of quantities from the set {A, B}. In this way we can 
learn, through the MC simulations, the partial derivatives that appear in the eigenvalue 
equation ([4.1|). Consequently, calculating the above derivatives at the fixed point and 
solving the eigenvalue equation ( [4.1D we obtain 



-^v — r 



' {N^rA' - {Nb)*b'V 



+ {NA)*n,{NBrA,, (4.6) 



which means that Xi, has been expressed in terms of quantities that all are measurable 
through the Monte Carlo simulations. Accordingly, we can learn the value of the critical 
exponent u through the formula (|3.4| ). We have applied this technique for a sequence 
of of the 3d SG fractals and in Table ^ we present our findings for 2 < 6 < 40, for 
the extended chain (in the bulk phase) fixed point {Ae,Be) together with the related 
critical exponent ue- As one can see from Table ^ the fixed point values for Be, decrease 
much faster then the values for Ae, when b increases. This means that we may neglect 
parameter B^"^^ (compared with A^^^) for larger b, that is for b > 12. In order to 
estimate the influence of the parameter B^^'> on the values of the critical exponent u, 
we calculated, by the MCRG method, critical exponent z/ for 6 = 12 with Be = 0. We 
obtained the following result u = 0.5989 ± 0.0003, which deviate 0.03% from the value 
u = 0.5987 ± 0.0003 (see Table |^) found using the nonzero value of B^'^\ Concerning 
the analogous results for the collapse transitions {6 critical point) in the bulk phase, 
we have to point out that the corresponding point {A0, Be) cannot be located because 
the initial part of the applied technique (the "homing procedure" ) does not possess the 
needed convergence. 

We now apply the MCRG method to find the crossover critical exponent (which 
determines the number of adsorbed monomers) for the polymer phase, on the 3d SG 
fractals, described by the symmetric fixed point ( p.8|) . To this end, we have to solve the 
eigenvalue problem for the second part (that starts with the third equation) of the KG 
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transformations ( p.8| )-( pTT2D , which reduces to solving of the equation 




0. 



(4.7) 



Here the asterisk indicates that all derivatives should be taken at the symmetric fixed 
point {Ae, Be, Ae, Ae, Be)- The above equation gives, in general, three eigenvalues for 
each b, but in practice it turns out that only one of them (to be henceforth denoted 
by X(f,) is relevant. Knowing we can determine the critical exponent (p through the 
formula 

Hence, in an exact RG evaluation of (f) one needs to calculate partial derivatives of sums 
( p.l0| )-( p.l2| ), and thereby one should find the coefficients ai{NA, Nb, Na-^, Na2, Nb^), 
02 (-^A, Nb, AOii, Na2, Nb^), and 6i(A^a, Nb, Na^, Na2, NbJ by an exact enumeration of 
all possible SAWs for each particular b, which has been accomplished for fractals with 
6 < 4. However, for 6 > 5, as in the case of the bulk phase, the exact enumeration 
turns out to be a formidable task. We have circumvent this problem by applying the 
MCRG method. Within this method, the first step would be to locate the symmetric 
fixed point. Fortunately, because of the structure of the symmetric fixed, the results 
given in Table ^ provide the complete information for this fixed point for a sequence of 
fractals with 2 < b < 40. The next step in the MCRG method consists of finding A,^, 
without explicit calculation of the RG equation coefficients. 

To solve the eigenvalue problem ( [4.7|) , so as to learn A^, we need to find the requisite 
partial derivatives. These derivatives can be related to various averages of the numbers 
NA-i^, Na2, and Nb^, of different SAW parts (of the types Ai, A2 and Bi) within a 



SAW path. Indeed, we may apply the relation ( |4.5|) where here X and Y stands for 
any pair of quantities from the set {Ai, A2,Bi}. Therefore, to calculate the derivatives 
( [4.5| ) for X, y G {Ai, A2, Bi} at the symmetric fixed point, one needs the nine averages 

{{NA.rA,, {Na2)\, {Nb.)\, {Na,)\, {Na2)\, {Nb^Ya,^, {Na^Yb'^, {Na^Yb',. {Nb.Yb'). 
which are all measurable through MC simulations. Solving numerically the eigenvalue 
equation ( |4.7| ) we obtain A^, and, finally, using relation (^78|), we find the values of the 
critical exponent (in the extended polymer phase), which are presented in Table 0, 
and discussed in the following paragraph. 

First, we would like to compare the results obtained, via the exact RG approach 
and through the MCRG technique, for the first three members (6 = 2,3,4) of the SG 
fractal families (given in Table |l| and Table ^ respectively). One can observe that 
the MCRG results for the critical exponents v and deviate at most 0.2% from the 
exact results. This very good agreement provides confidence in applying the MCRG 
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Table 2. The MCRG (2 < 6 < 40) results obtained in this work for the bulk fixed- 
point value parameters Ae and Be, and the SAW critical exponents ve and c/je for 
the 3d SG family of fractals. Each entry of the table has been obtained by performing 
10^ requisite Monte Carlo simulation. 

b {Ae,Be) ve 4>e 

2 (0.4311±0.0009, 0.0505±0.0023) 0.6742±0.0051 0.7484±0.0152 



3 


(0 


3421±0 


0004 









0021 





7148±0 


0037 


4 


(0 


2898±0 


0004, 


0.0122±0.0020) 





6414±0 


0012 





6901±0 


0036 


5 


(0 


2560±0 


0004, 


0.0067±0.0019) 





6315±0 


0010 





6707±0 


0028 


6 


(0 


2319±0 


0003, 


0.0038±0.0012) 





6239±0 


0009 





6496±0 


0024 


7 


(0 


2148±0 


0003, 


0.0020±0.0018) 





6169±0 


0006 





6333±0 


0019 


8 


(0 


2016±0 


0003, 


0.0012±0.0026) 





6130±0 


0005 





6198±0 


0016 


9 


(0 


1912±0 


0004, 


0.0007±0.0008) 





6087±0 


0006 





6023±0 


0017 


10 


(0 


1829±0 


0003, 


0.0005±0.0023) 





6048±0 


0003 





5894±0 


0013 


12 


(0 


1703±0 


0004, 


0.0001±0.0035) 





5987±0 


0003 





5686±0 


0012 


15 


(0 


1581±0 


0001, 


-) 





5933±0 


0002 





5385±0 


0010 


17 


(0 


1526±0 


0001, 


-) 





5899±0 


0002 





5213±0 


0010 


20 


(0 


1462±0 


0001, 


-) 





5869±0 


0002 





5014±0 


0012 


25 


(0 


1399±0 


0001, 


-) 





5817±0 


0002 





4666±0 


0008 


30 


(0 


1353±0 


0001, 


-) 





5804±0 


0002 





4325±0 


0008 


35 


(0 


1327±0 


0001, 


-) 





5759±0 


0001 





4275±0 


0007 


40 


(0 


1305±0 


0001, 


-) 





5755±0 


0001 





3765±0 


0007 



approach for a longer sequence of fractals (5 < 6 < 40). For the sake of a better 
assessment of the global behavior of the critical exponents u and 0, as a function of the 
fractal scaling parameter b, we depict the corresponding values (from Tables |l] and |^) 
in Fig. P and Fig. ^ respectively. One can see that u, being a monotonically decreasing 
function of b, crosses the narrow range determined by the predictions made for the three- 
dimensional Euclidean lattices, starting with u = 7/12 = 0.5833 [13|, passing through 
the values 0.5850 and 0.5874 |15|, and ending with u = 0.5882 |T6|. Crossing the 
interval of the possible Euclidean values for the exponent u occurs at the scaling fractal 
parameter b = 20. Interestingly, one may observe from Fig. ^ that the exponent 0, 
being monotonically decreasing function of b, crosses the estimated three-dimensional 
Euclidean value = 1/2 [|17| also at b = 20. 



5. Conclusion 



In this paper we studied the adsorption phenomenon of a linear polymer, in a good and 
bad solvent, on impenetrable boundaries of fractal containers modelled by the 3d SG 
family of fractals. Each member of the 3d SG fractal family has a fractal impenetrable 
2d adsorbing boundary (which is, in fact, 2d SG fractal surface) and can be labelled by 
an integer b {2 < b < oo). For the first three members {b = 2,3,4) of the 3d fractal 
family we have performed exact RG analysis. This analysis enabled us to establish the 
phase diagrams (for fractals with 6 = 3 and 6 = 4; the 6 = 2 case was studied previously 
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" 1/b 



Figure 5. Data for the end-to-end critical exponent v for the 3d SG family of fractals. 

The exact RG results arc represented by open triangles, while the MCRG results 
are depicted by solid triangles. The shaded horizontal band represents the region of 
estimated values (found in the literature) for the three-dimensional Euclidean critical 
exponent u. The error bars related to the MGRG data are not depicted in the figure 
since in all cases they lie within the corresponding symbols. 




" 1/b " 



Figure 6. Data for the adsorption critical exponent for the 3d SG family of fractals. 
The exact RG results are represented by open triangles, while the MCRG results 
are depicted by solid triangles. The dashed horizontal line represents the putative 
universal Euclidean value 1/2 of crossover critical exponent (j>. The error bars related 
to the MCRG data are not depicted in the figure since in all cases they lie within the 
corresponding symbols. 
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which turned out to be very rich from the physical point of view. Indeed, the 
phase diagrams disclose six different polymer phases that merge together at a multi- 
critical point, whose nature for 6 = 3 and 6 = 4, as well as the shape of the critical line 
that separates adsorbed from desorbed phases, depend on the value of the parameter t 
(associated with the monomer energy in the layer adjacent to the adsorbing boundary), 
which was not the case for 6 = 2. By analyzing the obtained phase diagrams we may 
conclude that similar diagrams can be expected for 6 > 4, but finding their exact pictures 
presently is not feasible. 

By applying the exact and Monte Carlo renormalization group (MCRG) method, 
we calculated the critical exponents z/ (associated with the mean squared end-to-end 
distance of polymers) and (p (associated with the number of adsorbed monomers) for 
a sequence of fractals with 2 < 6 < 4 (exactly) and 2 < 6 < 40 (Monte Carlo). The 
reliability of the MCRG results is manifested by the fact that in the cases 6 = 2, 3, 4 the 
MCRG results for u and deviate at most 0.2% from the exact results. Unfortunately, 
it was possible to implement this powerful (MCRG) method only in the case of the 
extended polymer phase. 

We find that, in the region studied, both u and monotonically decrease with 
increasing b (that is, with increasing of the container fractal dimension df), and the 
very interesting fact that both functions, z/(6) and 0(&), cross the estimated Euclidean 
values at, approximately, the same value of the scaling parameter [b ^ 20). Here we 
would like to point out that in the case of the 2d SG fractals both exponents u and 
also cross the corresponding Euclidean values, but not at the same value of 6 (at 6 ~ 27 



for u W^, while at 6 ~ 6 for <; 

On the whole, our findings should be useful in making the corresponding 3d models 
of the polymer adsorption phenomena in porous media. Besides, our results may serve 
as beneficial in constructing theories of the polymer adsorption phenomena for the 
homogeneous 3d lattices, where so far, to the best of our knowledge, an exact approach 
has not been made. 
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Appendix A. RG equations for the surface parameters in the case of the 
three-dimensional 6 = 3 SG fractals 

In this Appendix we give the exact RG equations for the surface RG parameters Ai, A2 
and Bi in the case of the 3d 6 = 3 SG fractal. We have found that these equations have 
the following form: 

A[=Y.Fl{Ai^A2,B,B,)A\ (A.l) 

1=0 
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(A.2) 



i=0 



(A.3) 



i=0 



where 



Fl {A,, A2, B, Si) = A? + 3At + AI + 2A? , 

F^^{Ai, A2, B, Bi) = B^AiAl{UA2 + IT^A^A^ + 132^iSi + 352>l?Si) , 
F^^(Ai, ^, Si) = 5^(22^^ + 132^2^251 + 484A?A25? + 484AfSf) , 

= 2 Ai ^ 9 A^ A\ + 10 Ax" A\ + 10 Ax^ + 8 B^ + 16 Ai^ As B^ + 

20 Ai^ A2 Bx + 12 Ai^ Si^ + 6 Ai^ Si^ + 32 A^^ B^ + 5^(6 A\ + 
14 Ai + 20 A^ + 64 Ai Al Bi + 32 Ai^ 5^ + 49 ^i Si^ + 
304 Ai^ Al Bi^ + 320 Ai^ Si^) , 

= A2^4AiA2^9AM^ + 16AM^ + 14VA2 + 3^2 + 9A>12 + 

12 Ai^ A^ + S Ai^ A2 Bi + 36 ^1^ A2 Bi + 40 Ai^ A2 Si + 4 Al B^ + 

24 Al Al Bi + 40 Ai^ Al Bi + 12 B^^ + 30 ^i^ Si^ + 
24 Al Al Bi^ + 88 A^^ Al B^^ + 48 Ai A2 Bi^ + 112 Ai'^ A2 Bi^ + 
128 Ai^ Bi^ + B{A Al Al + 12 Ai^ Al + S Ai^ Al + 16 ^i^ Al + 
4 + 24 Al Ai + 40 Ai^ + 8 Ai^ A2 Bi + 32 Ai^ A2 Bi + 
8 Ai^ A2 Bi + 12 Al Bi + 48 Ai Al Bi + 120 Al Bi + 24 Ai^ + 
64 Ai^ Si^ + 64 Al A^ Si^ + 176 Ai^ Al Bi^ + 176 Ai A2 B^^ + 
352 Ai^ A2 Bi^ + 264 Ai^ Bi^) + 5^(6 A^ + 12 Ai^ Al + 
36 Ai^ A^ + 20 A^ + 44 Ai A^ + 24 Ai^ A2 Bi + 24 Ai^ A2 Bi + 
96 Ai^ A2 Bi + 64 Al Al Bi + 232 Ai^ Al Bi + 36 Ai'^ B^^ + 
36 Ai^ Bi^ + 168 Al Al Bi^ + 168 Ai^ Al Bi^ + 928 Ai^ A2 Bi^) , 
= 2 A^ + 8 Al Al + 29 Ai^ A^ + 42 Ai^ Al + 34 Ai^ A^ + 4 A^ + 10 Ai A^ + 
14 Ai^ A2^ + 2 A^ + 28 Ai^ A2 Bi + 72 Ai'^ A2 Bi + 96 A/^ A2 Bi + 
4 A^ fii + 24 Al Al Bi + AA Ai^ A^ fii + 4 A^ 5i + 24 Ai^ fii^ + 
30 Ai^ Bi^ + 100 Ai^ Bi^ + 26 Ai Al Bi^ + 116 Ai^ Al B^^ + 
16 A^ Si^ + 64 Al A2 Bi^ + 112 Ai^ A2 ^i^ + 124 Ai^ Bi^ + 
136 Al Bi^ + 5(4 A^ + 12 Ai Al + 28 A^ A^ + 40 A? A^ + 
40 A\ Al + 2AAi + 48 Ai A^ + 76 A^ A^ + 20 A^ + 24 A^ A2 Bi + 
112Af A2 Si + 88Af A2 Si + 32 A^ Si + 136Ai Al Si + 192A^A^ Si + 
48 A^ Si + 24 Aj Si^ + 24 Af Si^ + 136 Aj Si^ + 96 Ai Al Si^ + 
768 Aj Al Si^ + 160 A^ Si^ + 448 Ai A2 Si=^ + 352 Aj A2 Si^ + 
944 A^ Si^ + 1320A^ Si^) , 
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F^^^2Al + 8 Ai Aj + 26 Aj + 36 Ai^ Al + 32 A^^ + 8 + 26 A^ A\ + 
32 Ai^ Ai + 24 Ai^ A2 Bi + 72 Ai^ A2 Bi + 80 Ai^ A2 Bi + 
SAlBi + 64 Al Al Si + 116 Al B^ + 24 A^"^ B^^ + 
24 + 100 + 74 Al Al B^^ + 

236 Ai^ Al Bi^ + 112 Al A2 Bi^ + 336 Ai^ A2 Bi^ + 292 Ai^ Si^ , 

F\^^ B{2Al + A A1AI + 12 Ai^Al + 6 Ai^Al + 6 AlBi + 12 AiAlBi + 
32 A^ Bi + 24 ^2 Bi^ + 44 Ai=^ Si^) + B^{2A Al + 
88 Al Al + 88 Al + 168 Ai'^ + 32 + 84 + 
88 Al Al Bi + 448 Ai^ Al Bi + 88 Ai^ Al Bi + 88 A^ Bi + 
132 Al Al Bi + 640 Ai" A2 Bi^ + 920 Ai Al Bi^ + 320 Al Bi^) , 

Fl^ = /I2 + 4 + 10 Ai^ A2 + 10 Ai^ A2 + 14 Ai^ A2 + 8 Ai^ A2 + 2 Ai^ Bi + 
8 Ai^ Bi+A Ai^ Bi + 18 Ai^ Bi + 5(8 Al + 24 Ai Al + 
36 Ai^ A^ + 32 Ai^ Al + 12 A^ + 12 Ai A^ + 36 Ai Al Bi + 
76 Ai^ Al Bi + 76 Ai^ A^ fii + 8 A^ Bi + 68 Ai A\ Bi + 64 Ai^ A2 Bi^ + 
88 Ai^ A2 Bi^ + 48 Al Bi^ + 88 Ai^ Bi^ + 304 Ai Al Bi^) + 5^(6 A2 + 
12 Ai^ A2 + 36 Ai^ A2 + 12 A/ A2 + 24 Ai^ A2 + 36 A^ + 56 Ai Al + 
104 Ai^ Al + 96 Ai^ Al + 50 A^ + 68 Ai A^ + 12 Bi + 
24 Ai^ Si + 12 Ai^ Bi + 12 Al Bi + 88 Ai Al Bi + 160 Ai^ Al Bi + 
296 Ai^ Al Bi + 72 A^ Si + 252 Ai A^ Si + 168 Ai A2 Si^ + 
216 Ai^ A2 Si^ + 168 Ai^ A2 Si^ + 160 Al Bi^ + 328 Ai Al Bi^ + 
328 Ai=^ Si=^ + 132 Al Bi^ + 952 Ai Al Bi^) 

Fj^ = 2 A2 + 4 Al A2 + 12 Al' A2 + 16 Ai=^ A2 + 16 Ai^ A2 + 12 Ai^ A2 + 8 A^ + 
24 Al Al + 42 Al' Al + 26 Ai^ A^ + 4 Ai' Si + 8 Ai^ Si + 8 Ai^ Si + 
20Ai^Si + 4A^ Si + 36AiA^ Si + 68Ai'A^ Si + 92Ai^A^ Si + 
20 Al A2 Si' + 76 Al' A2 Si' + 112 Ai^ A2 Si' + 120 Ai=^ Si^ + 
S(4 A2 + 8 Al A2 + 24 Al' A2 + 32 Ai^ A2 + 32 Ai^ A2 + 
24 Ai^ A2 + 22 A^ + 64 Ai Al + 100 Ai' Al + 86 Ai^ Al + 14 A^ + 
32 Al A^ + 8 Al' Si + 16 Ai^ Si + 16 Ai^ Si + 40 Ai"^ Si + 
6 A^ Si + 60 Al Al Si + 236 Ai' Al Bi + 208 Ai^ Al Si + 28 A^ Si + 
28 Al Al Si + 64 Al A2 Si' + 88 Ai' A2 Si' + 304 Ai^ A2 Si' + 
320 Al A^ Si' + 132 Ai''^ Si'^ + 132 Al Bi^) , 

Fj^ = 2 A2 + 4 Al A2 + 12 Al' A2 + 16 Ai^ A2 + 16 Ai^ A2 + 12 Ai^ A2 + 10 A^ + 
30 Al Al + 50 Al' Al + 40 Ai^ A^ + 8 A^ + 12 Ai A^ + 4 Ai' Si + 
8 Ai'^ Si + 8 Ai^ Si + 20 Ai^ Si + 6 A^ Si + 40 Ai Al Si + 
94 Al' Al Si + 112 Ai^ Al Si + 20 A^ Si + 58 Ai A^ Si + 
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20 Ai A2 Bi^ + 68 Ai^ A2 + 108 A2 B^^ + 36 Al B^^ + 
92 Al Al Bi^ + 100 Ai^ Bi^ + 32 Aj Bi^ + 192 Ai Al Bi^ , 

Fl^ = S(3 Ai^ Al + % Ai^ + 6 A^^ Al B^ + 44 A^^ A2 B^^ + 66 A^^ B^^) + 
(44 Al + 44 AiM^ + 22 + 320 Al Bi + 66 Bi + 
120 Al A^ Bi + 312 Ai^ A2 Si^ + 312 Ai^ A2 Si^ + 160 Ai A^ Si^ + 
968 Ai^ Si^ + 968 Ai Al Bi^ + 4308 Ai Bi^) , 

Fl^ = Ai^ A2 + 7 Ai^ A2 + 8 Ai^ A2 + 15 Ai^ Si + 8 Ai^ B^ + B{<6 Ai A^ + 

12 Ai^ Al + 12 Ai^ Al + 10 Al A^ + 16 Ai^ Al Bi + 60 Ai'^ Al Bi + 
24 Al At Bi + 44 Ai^ A2 Bi^ + 176 Ai^ A2 Bi^ + 88 Ai Al Bi^ + 
132 Ai^ + 576 Ai A2 5/) + B^{& Ai^ A2 + A/ A2 + 
54 Ai^ Al + 48 Ai^ A^ + 6 A^ + 28 Ai Al + 18 Ai^ Bi + 
48 Ai^ Bi + 128 Ai^ Al By + 136 Ai'^ Al Bi + 18 A^ fix + 
124 Al Al Bi + 432 Ai'^ A2 Bi^ + 318 Ai Al Bi^ + 132 Ai^ ^i^ + 
632 Ai^ Bi^ + 312 Ai Al Bi^ + 1936 Ai A2 Bi^ + 1452 Ai Bi^) , 
= 2 Ai^ A2 + 10 A/ A2 + 8 Ai^ A2 + 2 Al A^ + 14 Ai^ Al + 22 Ai^ Al + 
18 A/ 5i + 16 Ai^ Bi + 26 Ai^ Al Bi + 74 Ai^ Al Bi + 
12 Ai^ A2 Bi^ + 164 Ai^ A2 Bi^ + 44 Ai^ ^i^ + 120 Ai''^ B^^ + 
S(4 A1M2 + 20 Ai^ A2 + 16 Ai^ A2 + 4 A^ + 12 Ai Al + 
33 Ai^ A^ + 24 Ai^ A^ + 8 A^ + 36 Ai^ B^ + 32 Ai^ B^ + 
72 Ai^ Al Bi + 174 Ai^ Al Bi + 18 A^ Bi + 48 Ai A^ Si + 
88 Ai^ A2 Bi^ + 352 Ai=^ A2 Bi^ + 66 Ai Al B^^ + 66 Ai^ B^^ + 
592 Ai^ Si^ + 312 Al Al B^^ + 1452 Ai B^^) , 
= 2 Ai=^ A2 + 10 Ai^ A2 + 8 Ai^ A2 + A^ + 6 Al A^ + 17 Ai' Al + 
16Ai^A^ + 3A^ + 8AiA^ + 18Ai^Si + 16Ai^Si + 
26 Ai^ Al Bi + 72 Ai^ A^ Si + 7 A^ Si + 32 Ai A^ Si + 
44 Ai^ A2 Si^ + 202 Ai^ A2 Si^ + 64 Ai Al Si^ + 
66 Ai^ Si=* + 230 Ai=^ Si^ + 100 Ai Al Bi^ + 
44OA1A2 Si^ + 472 Al Sl^ 

Appendix B. RG equations for the bulk parameters in the case of the 
three-dimensional 6 = 4 SG fractals 

The exact recursion relations for the bulk RG parameters A and S in the case of the 
3d 6 = 4 Sierpinski fractal have the following form 

A' = A^ + 12 A^ + 62 A^ + 220 A^ + 782 A^ + 2426 A^ + 12 A^ S + 128 A^ S + 

776 A^ S + 3416 A^ S + 13324 A^ S + 18 A^ S^ + 72 A^ S^ + 432 A^ S^ + 



Adsorption of self -interacting polymers on fractals 27 

2736 + 13294 + 56004 A^ B^ + 19 A^ B^ + 1456 A^ B^ + 
8704 A'^ B^ + 48256 A^ B^ + 213968 A^ B^ + 36 A"^ B^ + 688 A^ + 
6880 A^ B^ + 33264 B^ + 173936 A^ B^ + 816 A^ B^ + 781456 + 
4350864 + 2904 B^ + 31616 A^ B^ + 248608 A^ B^ + 

2047360 B^ + 9995376 + 70400 A^ B^ + 779584 A^ B^ + 
5688896 A' B^ + 13632 B^ + 10369 + 562192 A^ B^ + 
640 A^ B^ + 1792 A^ 5^32432 + 250672 A^ B^ + 1552992 A^ B^ + 
7712 + 98528 5^ + 24273344 + 65568 A^ + 
354336 A^ B^° + 2960032 + 11284384 + 44440800 A^ + 

139392 A^ 5" + 6933440 A fi" + 25364736 + 20945408 + 

259424 + 3408768 A^ B^^ + 7491712 A^ B^^ , 

= + 24 A^ + 268 + 1624 A^^ + 6544 + 20288 A^'^ + 50676 A^^ + 

103904 y4^-5 + 173050 A^^ + 225108 + 215392 A^^ + 134968 A^'^ + 

42514 + 328 A^ B + 2584 + 13940 A^^ B + 56768 A^"^ B + 

188356 A^^ B + 491496 B + 1000000 + 1539632 A^^ B + 

1701920 A^^ B + 1192152 A^^ B + 397392 5 + 36 5^ + 

832 5^ + 14176 5^ + 85720 A" fi^ + 368404 5^ + 

1199040 5^ + 2954904 A^* B^ + 5370136 A^^ B"^ + 6814424 A^^ B^ + 

5299288 A^'^ B^ + 1901008 B^ + 2112 A^ B^ + 17424 A^ B^ + 

89344 A^° B^ + 482272 A^^ B^ + 1965376 A^^ B^ + 5900928 A^^ B^ + 

12721232^1^^5^ + 18817920 B^ + 16698368 B^ + 

6662824 A^^ B^ + 616 A^ + 1584 + 6388 A^ B^ + 

76968 A^ B^ + 523696 B^ + 2577592 A" + 9113852 A^^ B^ + 

23214816 A^^ + 40127080 A^^ B^ + 41398192 A^^ B^ + 

18548660 A^^ B^ + 12912 A^ B^ + 62560 A^ B^ + 451600 A^ B^ + 

2759488 B^ + 11688304 + 34977824 + 

70124768 A^^ B^ + 82557440 A^^ B^ + 41200784 B^ + 

1584 + 45056 A' B^ + 296744 A^ B^ + 2361856 A^ B^ + 

12623808 B^ + 43429664 5^ + 102376032 A^'^ B^ + 

139837056 5^ + 78454776 A^^ B^ + 7744 A*^ B'^ + 133440 5^ + 

1858832 A^ B'^ + 10467744 A^ B'^ + 46337232 fi^ + 

126642624 A" fi^ + 201751952 fi^ + 126633376 A^'^ B^ + 

17232 A^ B^ + 54192 A^ B^ + 950784 B^ + 8336624 5^ + 

42390624 A^ B^ + 139673656 5^ + 253941232 A^^ B^ + 

172479256 5^ + 431776 A^ 5^ + 6405824 A' B^ + 32770448 A^ B^ + 
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127799456 + 273263456 A^^ + 204194352 A^^ B^ + 

471328 A^ + 4395872 A^ + 21993584 + 
107371568 A^B^° + 262855104 5^° + 199987864 5^° + 
1224960 A^B^^ + 17655616^*^5^^ + 76932128 A'' B^^ + 
220464528 A^ B^^ + 151443088 A^ B^^ + 3702272 B^^ + 
11642752 A^ B^^ + 66511552 A^ B^^ + 148524304 B^^ + 
99652120 A^ B^^ + 3748096 + 30934336 A^ B^^ + 

110065872 A^ B^^ + 26465392 A^ B^^ + 6830208 A^ B^^ + 
30509600 A^ B^^ + 49942336 A^ B^^ + 20614528 A^ B^^ + 5153632 B 
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